Nonlinear N-parameter spacetime perturbations: Gauge transformations 
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We introduce N-parameter perturbation theory as a new tool for the study of non-linear relativistic 
phenomena. The main ingredient in this formulation is the use of the Baker-Campbell-Hausdorff 
formula. The associated machinery allows us to prove the main results concerning the consistency 
of the scheme to any perturbative order. Gauge transformations and conditions for gauge invariance 
at any required order can then be derived from a generating exponential formula via a simple Taylor 
expansion. We outline the relation between our novel formulation and previous developments. 
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In the theory of spacetime perturbations 0, 0, S 0, 0| (see 0, Q, Q for an introduction) , one usually deals with 
a family of spacetime models which, in most cases, depends on a single parameter A: M\ — (A4,{T\}), where M. 
is a manifold that accounts for the topological and differential properties of spacetime, and {7a} is a set of fields 
' on representing its geometrical and physical content (this formulation does not depend on the gravitational 
field equations). The parameter A labels the elements of the family and gives an indication of the "size" of the 
C*~) • perturbations, regarded as deviations of M\ from a background model Mq. It can either be a formal parameter, as 
in cosmology 15. l& Flc|, in back-reaction problems (see e.g. or m the study of quasi- normal modes of stars 

and black holes |l3t Il4|. or it can have a specific physical meaning, as in the study of black hole mergers via the 
' close limit approximation, in the analysis of quasi-normal mode excitation by a physical source, or in the modelling 
of perturbations generated by the collapse of a rotating star (see EM E3 an d references therein). 

There are however physical applications in which it may be convenient to use a perturbation formalism based 
£C) ] on two 5] or more parameters. For instance, in order to study general time-dependent perturbations of stationary 
axisymmetric rotating stars using a spherical background. In this case one can separately consider the stationary 
' axisymmetric rotational perturbations, for example up to second order in a parameter f2, then considering the coupling 
^j*. of these with the first-order time dependent ones (see 5] for further discussion) . As it should be clear from this example, 
the advantage of an N-parameter Non-Linear Perturbation Theory (NLPT) is that it allows us to make distinctions 
between different types of perturbations corresponding to different parameters, so that we can study their coupling 
and some non-linear effects without having to compute the whole set of higher-order perturbations. Such a framework 
may provide flexibility by allowing us to look at a given problem from different points of view. It may also allow 
us to choose a simpler (typically more symmetrical) background to model a given physical scenario. Given that, 
even in NLPT, the differential operators are those defined on the chosen background, this can drastically reduce the 
computations, and even change the way to perform them. 

The aim of the present article is to introduce a novel approach to the study of the gauge dependence of perturbations 
in NLPT which (i) deals with an arbitrary number of parameters; (ii) provides a closed formula for the action of a 
general gauge transformation, valid to any order; (iii) the construction and derivation of formulas of practical interest 
is simpler than in previous frameworks B, 0, . This new approach is mainly based on the application of the Baker- 
Campbell-Hausdorff (BCH) formula 19]. This has been used previously [20|. in the context of the back- reaction 
problem in cosmology, to derive second order one-parameter gauge transformations. Here we show how to make 
use of the full power of the BCH formula deriving the transformation between two given N-parameter gauges, each 
represented by an N-parameter group of diffeomorphisms, at an arbitrary order. Our formalism therefore also contains 
the conditions for gauge invariance for every perturbation order in N-parameters. 

We start by summarizing some relevant results regarding the mathematical structure of the single parameter NLPT 
(see 0, 0). The Taylor expansion of a general non-linear gauge transformation can be expressed in terms of Lie 
derivatives with respect to a set of vector fields which, order by order, constitute the generators of the transforma- 
tion. A closed formula for this expansion, valid at all orders, was found in Hj, using a new object dubbed Knight 
Diffeomorphism (KD), first defined in (cf. also [ill ISlfl- The analysis in Hani gives also the conditions for 
gau ge invariance at any given order, and provides the framework for the construction of gauge-invariant formalisms 
[23| . A similar construction has been attempted for the two-parameter case in [5|, where the action of a general 
gauge transformation on arbitrary tensor quantities was expressed in terms of the Lie derivatives with respect to a 
set of vector fields. However, since no natural generalizations of the KD approach were found, these expressions were 
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derived up to fourth order in the parameters by imposing, order by order, consistency conditions (see |24( for a related 
analysis and [25j for an application to gravitating strings). It must be pointed out that, although that derivation is 
not as elegant and compact as in or the one based on the BCH formula presented here, it still leads to the right 
formulas of practical interest, as we shall discuss. 

The basic assumption for the construction of a multi-parameter relativistic NLPT is the existence of a multi- 
parameter family of spacetime models = (Ai, {T^}), where Ai denotes the spacetime manifold and {T^} is a set 
of fields on Ai, describing their geometrical and physical content, which we assume to be analytic. These spacetime 
models are labeled by a set of N parameters A = (Ai, . . . , A N ) that control the strength of the perturbations with 
respect to the background spacetime model Mg, which describes an idealized situation. In order to construct the 
physical spacetime model as a deviation from the background model Afg, we need to establish a correspondence 
between them; what in the context of relativistic perturbation theory is called a gauge choice. This correspondence 
is established, for all A, through the action of a diffeomorphism of Ai: (p? : Ai — > Ai. The set of diffeomorphisms 
Q[ip] = {ip^ | A <E IR N } is chosen in such a way they constitute an N-parameter group of diffeomorphisms of Ai: 



<p: Mxl k 
(P,A) 



M 



<f(p, A) 



(1) 



The identity element corresponds to A = 0, i.e. ip^(p) = p. Moreover, a consistent perturbation scheme should 
have the property that perturbing first with respect to a given parameter, say A P , and afterwards with respect to 
another parameter, say A Q , must be equivalent to the converse operation. We can implement this idea by imposing 
the following composition rule for the group Q[ip]: 



V A , A' , n o <p x , = p x 



(2) 



This property implies that the group is Abelian. It also implies that we can decompose ip^ into N one-parameter 
groups of diffeomorphisms (flows) that remain implicitly defined by the equalities (we have N! equalities) 



<P\ = <P(\ U 0,...,0) ¥>(0,A 2 ,...,0) ° • • • ° V(0,0,...,A N ) 

= V(o,A 2 ,...,o) ° f(\ u o,...,a) ° ' ' ' ° ¥>(o,o,.. .M) 



(3) 



The action of the flow 



<^ Q Am q ^ is generated by a vector field, C M (M= 1, . . . ,N), acting on the tangent space of 
Ai, and the Lie derivative of an arbitrary tensor field T with respect to C M is given by 



£f T = 



d 



T 



Am=0 



where the superscript * denotes the pull-back. Since the group is Abelian, the vector fields £ M commute 

[C P ,C Q ]=0 (P,Q = 1,...,N). 
The Taylor expansion around A = of the pull-back associated with the flow ip^ Am ^ is given by 



V(0,...,A M ,...,0) 



oo . i. 

A* 



fc=0 



This expression can be written in a more compact way by using the formal exponential notation 



V>(o,...,x M ,...,o)T = exp(A M £ CM )T = e A ^<M T , 



(4) 



(5) 



(6) 



(7) 



which provides a clear operational way for working with groups of diffeomorphisms. From expressions © and we 



can derive the Taylor expansion of the pull-back ip^T 



fci,...,fc N =0 \p=1 p ' / 

Using the exponential notation and the commutation relations (0 we can write it as follows: 



(8) 



JJexp(Ap£ Cp ) 



T = exp ( Ap£(- 



T . 



(9) 
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In a given gauge <p, the perturbation of an arbitrary tensorial quantity T is defined as 

AT X : =^ T X- T o- (10) 
The first term on the right-hand side of (|10|) can be Taylor-expanded around A = using © to get 

fci,...,fc N =0 \p=l 



where fc = (fci, . . . , fc N ) and 



^felH hfeN 



fei a\fcN *A J 



d\f ■■■ax 



U £ Z t > ( i2 ) 



p=i 



which defines the perturbation of order (k\, . . . ,k N ) of T (S^T := Tg). The total order of a perturbation can be defined 

as n,T '■= k\ + ■ ■ ■ + fc N . As a consequence of these definitions we have that AT-f and S^T are fields that belong to the 
background spacetime model Mg. 

Let us consider two different gauges <p and -0, with generators . . . , V ( N ) and . . . , ^C N ) respectively. For all 

A, the objects defined in these two gauges can be related by a diffeomorphism : M — ► M given by 

A 

$ A : =^A l0 ^A=^-A 0?/, X- ( 13 ) 

This is what is called a gauge transformation in perturbation theory. The family of all the possible gauge transfor- 
mations for two given gauges if and ip 

$ : M x TR N — ► M 

(p,X) | — > $(p,A) = $ x (p), (14) 

is not in general a group of diffeomorphisms. The action of the gauge transformation <£>^ can be written as 

= exp f^A P £« fp J exp f-^A P £<^ p jT. (15) 

Using group theory techniques we can write (|15[) as the action of a single exponential operator. This can be explicitly 
done by using the BCH formula [l0j (cf. 20] ). This formula, which can be applied to any two linear operators A and 

B, reads 

e A e B = e f(A,B) 5 (16) 

ft + qi > 1 

where the following notation has been used: 

[Xi* 2 *3 • • • X n ] :=[... [LY 1; X 2 ], X 3 ], • • • , X n ] . (18) 

Then, the BCH formula can be seen as an expansion in commutators of the initial operators A and B. Up to two 
commutators this expansion is given by 

/(A B) = A + B + \ [A, B] + ± [[A, B], B] + A [[B, A], A] + ■ ■ ■ . (19) 

This infinite expansion can be truncated and becomes finite when some commutators vanish (for a solvable Lie 
algebra). For instance, if [[A, B],A] = [[A, B], B] = 0, then the BCH formula only contains the three first terms shown 

in on. 
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The application of the BCH formula, (|16fl and (|17fl . to the construction of the gauge transformation <f>^ i|15|) 
constitutes the main point in our approach. As we are going to see, it provides an operational apparatus to compute 
all the perturbation orders as well as expressions for the vector fields that generate the transformation between different 
gauges. This supposes an important advantage with respect to the approach to two-parameter NLPT considered in 
[5j, which is based on a construction order by order and no close expressions to every order can be obtained. In what 
follows, we show how to use the BCH formula to obtain closed expressions at every order, in particular for the gauge 
transformation generators. 

To apply the BCH formula to Eq. (|15|l we have to choose A and B in (|16|) and l|17|) as follows: 

N N 

A = ^A P A Cp and S = -^A P £^ Cp . (20) 
p=i p=i 

Since A and B are linear in the parameters, the number of commutators in a given term in the expansion of f{A, B) 
coincides with the total perturbation order n^. 

Using the properties of Lie derivatives we can then express the gauge transformation $^ in the following way: 



(21) 



This can be rewritten as: 




<t>* x T = exp{ ^ |n?r|/( e -lh (22) 



where £^ s denotes the identity operator / and the rest of terms are Lie derivatives. This is a consequence of the fact 
that A and B are linear combinations of Lie derivatives (|20|l and that the functional f(A, B) is a linear combination 
of A, B, and commutators formed out of A and B H19(l . Then, this introduces an infinite set of vector fields 
{£j | k 6 IN N — {0} }. By direct comparison of 1(19120(1 and (1221 we can find the explicit expressions of these vector 
fields £2 directly in terms of the generators of the gauges <p and ip. 

We can then derive an expression for the gauge transformation up to a given order in the perturbation parameters 
by simply expanding the exponential (|22|l . Up to third total order (tit — 3) we obtain 

N N , s 

$tr = t + ^a p £ s . t + - a p a q [ £ e- - + £ e- £ e- 1 T 

P— 1 P,Q— 1 \ f w r 



r J2 a p a q a r [ £ e s 



3 

0/ v ' "P ' "Q " "K \ " t - - - 1 o 
. V k T> + k n + k T, 2 



P,Q,R=1 



pT~Q T ~ R a ^ fc P + fc Q fc R 



% % % T + 4 (A), (23) 



p-i 



where k P = (0, . . . , 0, 1, 0, . . . , 0) and {A, B} denotes the anticommutator of A and B. From (|23(l we can easily derive 

4 T - 6 P 



the transformation of a given perturbation from the gauge ip to the gauge ip, 8\T — Sj\T. First, from 1(13(1. the 



pull-backs ^iT and V^T 1 are related by 

^ T X = $ M T X' ( 24 ) 

Then, using 1)10(1 and ((11(1 we have that 

00 / N A fep \ - 
fei,...,fc N =0 \p=1 P / 

and for ip^T^ we only have to replace ip by if). From l(23|l - 125(1 we can extract the expressions for the f^T's in terms 
of Lie derivatives of the c^T's. In the particular case N=2 and k = (1, 1) we find 

,(M) T = S £A) T+££n xm T+£ein x£,o) T 



To . (26) 
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With the aim of comparing the formulation here introduced with previous approaches to NLPT, we will show now 
how to recover standard one-parameter NLPT. Let us consider two gauge choices: ip and ip. For a given A, the action 
of their associated pull-backs can be written in the exponential notation as: (p* x T — e x£v <T, and ip^T = e *<T. A 
gauge transformation between these two gauges is then described by $a = 'Px 1 Using the exponential notation, 

its action can expressed as follows 

$ * r = ^£i, te -XU iT (27) 
The result of using the BCH formula can be written as 

( °° A™ \ 

*i T = ^\Eri £ e~) T ' (28) 

where {£„ | n € IN — {0}} is a set of generators of These can be expressed in terms of the gauge generators % and 
. From l|19|l the three first generators are: 

& = *c - *c, & = [x % , 6 - \n + % rc, *ci] . (29) 

Up to third order, (|28|l gives [the N=l subcase of : 
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3! 



(£ 6s + I , £ i2 } + 4) T + 4 (A) . (30) 

This form of &* X T , derived through the BCH approach, is not the same as the one obtained using KDs in 0)0]. 
However, as we are now going to show, the resulting gauge transformations are - order by order - equivalent, as 
expected, since both cases are expansions of the same exact expression Q27|l . The KD is defined as 0, E| : 

where the (p^ are one-parameter groups of diffeomorphisms (flows). The main idea behind this concept is that a 
family of diffeomorphisms {<&\ | A G ]R} can be approximated at a given order in A, say k, by a KD of order k. 

Therefore one can approximate $a by in the following sense Q: 

$ » T _ $ (fc)* T = Qk+l ^ _ ^ 

(k) 

The action of the pull-back of $^ can be expressed, using the exponential notation, as 

$W*T = e Air *i ■■■e 2kr£ *kT, (33) 

where {Xn} n= i k 1S ^ ne se ^ °f generators of the family <£>( fc ), and each Xn is the generator of the flow tjy- n ' . Like the 
£n's, they can be expressed in terms of the gauge generators v ( and Hence, we can find the relations between the 
£„'s and Xn's. Up to third order we have 

3 

6=Xi, 6=X2, & =X3+ gbciiXa]- (34) 
Therefore, the expansion for that we obtain from the expansion of $^*T is (up to third order): 



A^ 
3! 



(£ X3 +3£ Xi £ X2 + £l i )T + OHX), (35) 



i.e. the result in P|. |_|. Comparing the expansions (|30|l and (|35|) we see that they have different structures. However, 
substituting Xli Xi an< 4 X3 from (|34|l into l|35|) we obtain (|30|) : thus, these two expansions of the gauge transformation 
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(|27|l are equivalent. Our formulation, Eq. (|30|) . leads to an expansion with terms of the form • ■ • £g • • ■ • • • T 
with I < k , which do not appear in the formulations of [due to the ordering introduced by the KDs, see i|31[l ] 
and Comparing further our formulation with the order by order approach in Q we see how the use of the BCH 
formula naturally selects a unique expansion for the gauge transformation, in contrast with ||, which contains freely 
specifiable parameters. 

To sum up, we have presented a formulation of N-paramctcr NLPT in which we have a (unique) closed formula 
for the expansion of general gauge transformations and whose consistency is given by construction, shedding light 
onto the underlying mathematical structure. The importance of this result is even more clear if considered in the 
context of practical applications of relativistic perturbation theory, where the issue of comparing results obtained in 
different gauges and the related problem of constructing gauge-invariant formulations have always been crucial to 
obtain physically tr ansp arent results Q, |2^ | . These issues become even more important when dealing with non-linear 
perturbations jj, 0, l22| and more than one parameter. Our formalism provides the gauge transformations and the 
conditions for gauge invariance for every perturbation order in N-parameters (explicit conditions for the 2-parameter 
case where given in 0). 

In retrospect, one may wonder why our general results have not been previously derived, given that the BCH 
formula has long been known in differential geometry |l9j . The likely answer is that, although the problem of gauge 
dependence is as old as relativistic perturbation itself untill recently spacetime perturbations have mostly been 
considered at first order only and for a single parameter, and consequently gauge transformations have always been 
dealt with at the most elementary linear level, where the BCH formalism and the exponentiation Q on which it is 
based are superfluous. When the problem of gauge transformations has been considered in NLPT for the case of one 
parameter, two routes have been followed. In :S, 4] KDs have been introduced and used (see also [llj for an equivalent 
second order result and |2l| for some basic fomulas), and in particular in Q a closed formula was derived to generate 
gauge transformations at arbitrary order. In [U on the other hand the BCH formula was used, for one parameter 
at second order. As we have illustrated above in the one parameter case and up to third order, the two routes are 
equivalent in that they provide equivalent gauge transformations at the required order. On the other hand, the gauge 
transformations derived in [f| contain freely specifiable constants (linked by sets of constraints) that are not present in 
the BCH derivation presented here. Again, order by order the gauge transformations are equivalent, with one specific 
choice of the constants corresponding directly to the BCH derivation, and other choices connected by appropriate 
relations between the two sets of generators of the gauge transformations. From the point of view of generality 
elegance and compactness of the derivation the BCH approach presented here is by far superior to that followed in 
@. However, for practical purposes one is interested in the gauge transformations at a given order, e.g. (|2*S|l . and in 
this perspective we believe that the formulas with freely specifiable constants in [j| may still be useful. Indeed, the 
typical problem (see e.g. 0,0) * s that one wants to know how to transform between two pre-assigned gauges. In this 
case the unknowns of the problem are the generators of the transformation. Then, one faces integration calculations, 
and given that two different choices of the freely specifiable constants correspond to integration in a different order, 
it may turn out that one specific choice of constants is better in solving the problem. 

We want to finish by discussing the potential applications of N-parameter NLPT. First of all, it is important to 
remark that perturbation theory in general relativity, and in other spacetime theories (some of them of great relevance 
nowadays), remains the main alternative to fully numerical methods in a context in which one has to deal with sets of 
non-linear field equations. Depending on the physical problem at hand, it is sometimes necessary to go beyond simple 
linear perturbations, considering higher-order contributions. In this sense, a multi-parameter perturbative scheme as 
the one presented here allows us to select only the higher-order perturbative sectors relevant for the physical problem 
under consideration, simplifying in this way the calculations involved. 

There are already quite few applications of the one-parameter NLPT at second order in the literature. In cosmology, 
the evolution of perturbations in two different gauges is explicitly worked out 0, E| an d applications to the Cosmic 
Microwave Background have been considered 26]. Further applications in the cosmological context can be found in 
[Tol l2(il ] . In an astrophysical context there are a number of studies of sources of gravitational radiation: in 0, |53 
oscillations during gravitational collapse have been analyzed; in |15| | the so called close limit approximation is used 
to study the outcome of black hole mergers; in [2^ a second-order gauge-invariant perturbative scheme for the Kerr 
metric has been developed. The N-parameter NLPT opens the door for new applications in spacetime theories (see 
|25) for an example). In the general relativistic case, it can provide a new way of studying slowly rotating relativistic 
stars, and it can be the main tool to study other issues as for instance the non-linear coupling of oscillation modes 
of relativistic stars, or in cosmology to study the combined effect of magnetic fields and linear perturbations in the 
matter distribution. 
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